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THE  ETOT  OP  LENGTH  OH  THE  AERODTNAKIC  CHARACTERISTICS  Of 
BOOTS':  OP  REVOLUTION  IN  SUPRSONIC  FLIGHT 


ABSTRACT 

As  a result  of  a joint  interseiTice  rasa  arch  effort,  the  a ray 
Easy  Spinner  Heckat  program  eonaiating  of  asn  than  three  hundred 
■odels  of  serious  lengths  has  bean  fired  on  Btf.'a  precision  Ttm 
Plight  Spark  Range.  The  data  obtained  froai  these  firings  are  analy- 
sed to  pnvrlde  a good  datendnatlan  of  the  effot;’  of  aodel  length 
on  the  aerodpasado  coefficients  for  supersonic  Mach  Barbara.  The 
effect  of  length  an  djrranta  stability  la  c-nsidsred  in  detail. 

Appendices  provide  a auanarr  of  theoretical  relations,  contention 
relations  between  aarodptuodo  and  ballistic  nceeaclatura,  and  a loll 
tabulation  of  the  aaperlaantal  data. 
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INTRODUCTION 

In  19li6  «n  intensive  program  investigating  Lhe  dynamic  properties 
of  spinning  missiles  1a  parsonic  flight  uaa  instituted  aa  a Joint 
interservioe  affort.  The  Arey  throu^di  the  Free  flight  Aerodynamic  a 
Pranch  of  the  Ballistic  Rc«sirch  Labors tcrlf  s a preed  to  aeuuf*cture, 
measure,  and  fin  tla  one  hu-dred  and  t urtr-f3*»  r squired  by 

the  program  .rnd  thr  Nary  agreed  to  contribute  the  plats  measurement 
and  data  redaction  facilities  of  its  ccmputing  group  under  Dr.  Z Isaak 
Kopal  at  the  Hgmeetuswtts  Institute  of  Technology.  A- report  on  the 
dat-s  reduction  process  was  issued  by  the  M.I.T.  poop. 

The  basic  Interest  of  ths  prograa  was  la  ths  dynarrle  stability 
of  spinning  bodies  of  ravolutd rm  with  an  erphasis  oa  configurations 
possessing  large  On areas  ratios.  Of  the  six  ballistic  coefficients 
of  the  Kelly-McShane  linearised  theory  tdxich  affect  dynamic  atafcj 
rary  few  aeasurenenta  of  the  Karats  acnent  coefficient,  ”.( 


the 


dating  aanent  coefficient,  1^(0^  ♦ ),  had  boon  aada 


V 


Information  on  these  coefficients  was  therefore  to  be  an  important  result 
of  thw  pr-v.va  Since  these  quantities  are  functions  of  the  cen-er  of 
mss  location,  th«»e  force  coefficients  which  characterise  these  Amotions 
also  had  tc  bn  found.  Thaaa  were  the  Magma  fnron  coefficient. , 

VC»J  , the  detqring  force  coefficient,  ),  end  the  normal 

pc  q d 

force  coefficient,  £.(CV  ).  The  rsaainlng  three  ballistic  coefficients 

“• 

are  the  overturning  iaooent,  K^(Cy  ),  which  ia  of  pries  imports:  .e  to 

4 

gyroscopic  stability  and  Is  essential  for  consideration  of  dynamic 
stability,  and  the  two  axial  coefflciauter  axial  drag,  X^fC^),  and 

spin  deceleration,  E^(C^). 

the  determination  of  most  of  these  ballistic  coefficients  especially 
h is  usually  quite  difficult  it,.'  wind  tinnsls  and  hence  the 

Selection  of  mb's  Tree  flight  opart:  Range  for  the  a tuny  was  a logical 
choice.  This  range  «t  precant  consists  of  forty-el r enark  etat-lonu  which 


1.  Xopal,  Kavartagh,  and  Roller,  A Manual  of  Reduction  of  Spinner  RooIgu 
Saadcw grama.  Center  of  Anaiys^-*,  Iec£.  Rept.  JsoTTT  vout  uf  print}. 

2.  the  symbols  appearing  in  parenthesis  after  the  ballistic  c.wfricient 
are  the  corresponding  aerodynamic  coefficient.  See  Appendix  » for  a 
discussion  of  the  precise  correspondence.  Hrackets  will  be  used  to 
identify  the  uumber  of  individual  publications  listed  in  the  References. 
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are  distributed  over  a distance  of  two and  nighty  feet. 

Figvr-  1 shows  a view  looking  down  rang*' with  the  spa-k  cylinders 
on  the  left  and  the  plate  holders  on  the  right.  The  brass  squares, 
located  at  each  station,  shield  single  loops  of  wire  which  font  part 
of  the  electrostatic  tr-.ggering  circuit.  As  the  gun  launched  models 
peas  each  of  the  spark  stations  accurately  located  to  0.001  feet, 
the;  are  photographed  simultaneously  la  the  vertical  end  horizontal 
plants  by  a short  duration  spark  discharge.  The  time  of  occurrence 
of  the  discharge  for  ten  o l hue  stations  Is  measured  to  an  accuracy 

of  10”°  seconds,  from  the  photographs  the  spatial  coordinates  of  the 
alas  He  are  obtained  to  an  accuracy  of  .001  feet  in  position  and 
three  minutes  of  are  in  angular  orientation.  The  ballistic  coefficients 
ere  computed  from  these  data.  QlOj  • 

When  the  program  was  originally  set  up  it  was  decided  that  five 
rounds  were  necessary  to  determine  one  value.  Since  it  was  planned  to 
obtain  reliable  values  of  the  force  coefficients  from  the  yawing  motion 
of  identical  shapes  possessing  different  centers  of  mass,  this  resulted 
in  the  requirement  that  three  different  canter  of  mass  positions  for 
each  shape  be  fired.  In  order  to  stu^y  stability  over  a reasonable 
range  of  supersonic  Mach  numbere  data  were  to  be  obtained  at  three 

Mach  number e , 1.3,  1..9,  and  2.$.*  These  considerations  meant  that 
forty-five  rouoo  i would  be  needed  to  complete  the  study  of  each  design. 

Body  length  was  selected  to  be  the  basic  design  variable  and  the 
shape  was  to  b«  representative  of  service  spinner  rocket  dacigna.  Tar 
these  reasons  the  canon  head  shape  was  choson  to  bs  a secant  ogival 

head  taro  calibers^  long  whose  radius  of  ogive  was  twice  the  tangent 
ogival  radius.  Three  body  lengths  of  three,  five,  and  seven  calibers 
were  agreed  upon  and  thus  find  the  program  sise  at  135  mode,.".  a 
drawing  of  the  configuration  with  the  actual  center  of  mass  location 
indicated  is  shown  in  Figure  2 end  e tabulation  of  the  physical  character- 
istics of  the  Models  la  given  In  Table  Cl  in  Appendix  C. 

In  order  to  be  able  to  make  flow  computations  e smooth  contour  was 
required  and  hence  no  rotating  band  was  used.  Spin  was  imparted  by 
means  of  a p re -engraved  aluminum  sabot  placed  immediately  behind  the 
model.  At  first  friction  coupling  was  employed  but  later  it  was  found 
to  be  necessary  to  connect  mcdel  and  sabot  by  means  of  a cruciform  key. 

1.  See  [8}  and  £9}  for  descriptions  of  the  range.  Appendix  A offlcQ 
giv-ss  i more  recent  d:  scrip  tier.  *»f  ‘■.he  -an  go. 

2.  At  the  time  the  intermediate  Mi  h nutr.oer  firings  were  made,  no  suit- 
able gun  vaa  available  and  hence  the  filings  were  made  at  the  sligitly 
lower  Mach  number  of  1.8.  Later  using  a gun  contributed  by  the  Navy, 
it  was  possible  to  lire  at  the  high  Mach  number  of  2.1?. 

3.  A caliber  is  one  diameter.  In  this  program  20r » models  were  used  snd 
hence  one  caliber  correiipondsd  to  20nm. 
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Finally  a copper  obturating  sup  was  placed  behind  the  sabot.  A nodal 
tilth  its  key,  sabot,  and  copper  cup  is  shown  in  Figure  3.  Althougi 
these  av.  Till  ary  launcuini*  c opponents  usually  separated  quite  quickly 
from  the  model  due  to  their  1 igh  drag,  for  some  rounds  it  wao  possible 
to  observe  this  separation.  (Soe  Fig.  4).  In  hi  tbs  missile  is  about 
two  feot  from  the  gun  and  none  of  the  corponents  dare  separated.  The 
frtgnects  »bo7e  the  miweils  are  Arcs  paper  wadding  used  in  the  cartridge 
case.  Figure  h-s  diows  the  cottier  cup  s-parating  at  fire  feet  fron  the 
gun.  At  twenty-three  feet  (lC)  the  sabot  is  clear  of  model  and  the  key 
is  dropping  free.  Finally  the  sabot  separates  to  over  four  calibers 
and  the  •ieeile  is  in  free  flight  twenty-eight  feet  frost  the  gun,  (ud)« 

In  the  firing  portion  of  'he  program  considerable  difficulty  was 
experienced  !*>  launching  the  models  properly.  In  order  to  obtain  the 
aerodynamic  moment  coefficients  a yawing  motion  whose  magnitude  is 
between  half  a degree  and  six  degrees  is  required.  The  lower  limit  is 
imposed  by  measuring  accuracy  while  the  upper  indicates  the  limitations 
of  the  linearised  heory.  ( 

In  order  to  obtain  the  program  requirements,  330  models  were  ' 
launched.  It  was  possible  to  perform  a complete  yaw  reduction  on  109 
of  these  of  vhicn  V nao  average  squared  yaw  of  over  30  and  so  were  set 
Included  In  tne  rnalysis.  The  overturning  moment  coefficient  was  ob- 
tained fr-.a  19  additional  rounds  whose  yawing  motions  were  too  email  for 
a collets  reduction.  Of  the  models  launched  145  were  suitable  for  drag 
reductions.  A total  of  162  rounds  provided  .usable  data.  A good  indi- 
cation of  the  improvement  In  firing  efficiency  as  experience  was  gained 
is  shown  by  Figure  5- 

The  body  of  this  report,  which  is  basically  concerned  with  the 
data  resulting  f*um  these  firings,  will  be  divided  into  three  aajo-.* 
parts.  The  first  part  will  deal  with  axial  force  and  moment  coeffi- 
cients, the  sacond  with  transverse  ft-ce  ^ad  moment  coefficients,  and 
the  third  will  analyse  the  dynamic  stability  characteristics  of  the 
configurations.  In  the  appendices  all  theoretical  relations  are  stetea 
end  the  conversion  from  the  ballistic  nomenclature  to  aerodynamic  is 
derived. 

AXIAL  FOflCB  AND  MOMENT  COEFFICIENTS 
(a)  Drag  Fcrco  Coefficient 

Denoting  tie  component  of  the  aerodynamic  force  along  the  a:d.s  of 
the  missile  fcy  F^,  the  udal  drag  coefficient  can  then  tv  defined  by 

the  equation: 

T\ 

(1)  E_.  i- 


~DA 


P 


vhoro  p » air  density 
d « diameter 
u^»  axial  velocity 
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If  F^  in  (1)  is  replaced  by  Fy,  the  component  of  the  aerodynamic  force 
directed  along  the  trajectory*  the  more  familiar  drrg  coefficient 
Kq  ■ tj  Cy  is  defined.  If  cross  spin,  p,  „»  neglected,  it  can  be  shewn 
that  Kjj  - Xjjj  cos  C ♦ K^S  sin  6 ■ j 62  where  6 is 

the  magnitude  of  yaw  and  Kg  is  the  normal  force  coefficient.  The  dra^ 

coefficient*  along  with  the  Mach  number,  is  determined  for  each 
round  by  means  of  a polynomial  least  squares  fit  of  time-position 

data.  filjJ 

Since  Sjj  is  a function  of  several  parameters,  basically  shape  cf 

missile,  Mach  number,  M,  and  the  squared  magnitude  of  the  yaw,  the 
effect  of  the  different  parameters  have  to  be  separated.2  The  assmption, 
therefore,  is  made  that  the  drag  dependence  on  yaw  is  linear  in  the 
average  squared  yaw,  6 , for  the  particular  round.  This  average  is 
over  the  distance  between  the  first  and  last  timing  stations  and  can 
easily  be  confuted  from  the  parameters  of  yaw  reductions.  Drag  values 
were  grouped  in  intervals. of  Mach  number  Iocs  than  0.1  long  «ud  ths 
parameters  :f  th's  linear  function  for  each  group  found  from  a least 
squares  ft  of  the  drag  coefficients  and  corresponding  mean  squared 
yaws.  The  zero  yaw  values  of  the  drag  coefficients,  Xq  , were  thoa 

ccqputed  and  use  via  made  of  the  empirical  Q function  derived  by  - 
Thcaes  , JJl2J  , to  derive  their  dependence  on  Mach  number. 

The  Q function  essential!-'  assumes  an  inverse  quadratic  dependence  of 
Kp  on  Mach  number  and  can  be  written  in  the  form 


- a ♦ bM 


where  a and  b are  empirical  constants.  This  function  provides  a 
good  description  of  the  drag  coefficient's  dependence  on  Mach  nuabar 
for  supersonic  Mach  numbers.  The  constants  a and  b for  each  model 
length  were  determined  by  the  usual  least  squares  method. 

Since  the  yaw  drag  coefficient  Kp  was  contaminated  by  Mach  number 

S 

effects,  the  rounds  of  each  interval  were  converted  to  central  Mach 
numbers  by  means  of  (2)  and  new  yaw  drag  oeflicients  determined. 

2ne  process  cm  noy  )>**  iterated  and  convergence  is  rapid.  Table  I 


1.  Throughout  the  report  we  will  indicate  for  each  ballistic  t its 
exact  relationship  to  tie  corresponding  aerodynamic  C’s  (See 
Appendix  B for  details  of  the  conversion,  process.) 

2.  in  '&■  dspendei.ee  of  on  spin  which  has  been  measured  it 

described . . 

3 In  £l!?J  the  technique  for  handling  drag  data  which  ic  employed 
here  i3  described  and  applied  to  a study  of  the  effect  oh  drac  of 
systematically  blunting  the  headshape. 

3 
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preseu.s  the  fintl  values  of  the  a's  and  b*s  together  •.-•ith  the  Kp  2 
used.  (These  *M2  may  he  consorted  from  1 /sq.  radian  to  l/sq.  de|r».e 
by  the  factor  3280  j3~g|~£  ) Dse  defining  relation  for  Kg  2 is 


(j) 


fD  * \2  «* 
o 6 


TABLS  I1 


1.  All  errors  calculated  in  this  report  are  standard  errors.  In  order 
to  consort  to  probable  error,  the  mil  tip  lica  tire  factor  0.6/U5 
should  be  employed. 

2.  Si  ice  Kjj  is  about  1.0  this  table  indicates  that,  the  axial  drag 

coefficient,  K-.  i z also  otr  :ngV  dependent  on  tho  magnitude  of 
yaw.  ** 
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1 1 

u 


; 

- lAuaction  Parameters 

f 

5 cal. 

7 cal. 

9 cal. 

y 

.921* 

a 

.927 

.:*> 

b 

.157 

.161* 

.172 

V 

1 2 
Square  radian? 

I-' 

Mach  No 

5 cal. 

7 cal. 

9 cal. 

1.3 

3.1  ♦.& 

1*.7  ♦ .6 

3.1*  ♦ .1* 

' ■<  1 

l'  i 

1.8 

2.6  „>3 

2.6  ♦ .2 

2.3  * .3 

r‘. 

r ■; 

2.5 

2.7  ♦!! 

2.7  ♦ .3 

2.9  +1.0 

\ ; 
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Tus  wlgxxal  X^'s  for  each  round  and  their  mean  >qtuu*ed  ft* a 

are  tabulated  In  Appendix  C.  An  indication  cf  theU'  internal 
statistical  *?curacy  is  provided  by  the  representative  standard 
errcre  provided  by  Table  C7«  These  standard  errore  of  about  O.SS 
seen  to  conflict  with  the  actual  dispersion  of  about  }%  as  seen  in 
Figures  9,  10,  and  11*  This  rp  .ead  la  quite  easily  explainablo  by 
tho  wiae  variati  n o*  boundary  layers.  Figures  6 and  7 clearly 
il Insert ic  this  points 

A Airthar  check  on  this  ixplans.tton  can  be  obtained  by  Beans 
of  a siuple  esaputatisn.  It  : » p ossi'  la  to  ccepute  approximate 
values  of  the  skin  friction  drag  coefficient  K^gy  for  both  turbu- 

lert  and  laminar  boundary  layers.^"  The  results  of  this  computation 
are  listed  in  Table  IT. 

TABLE  n 

/ 

Skin  Friction  Drag  Coefficient,  x 10 


H 5 cal*  7 cal*  9 c*l» 


turb. 

laa. 

% diff.2 * * 

turb. 

lam* 

* diff. 

turb. 

lam. 

% diff, 

1.3  2.1 

.5 

ID 

3.0 

Ha 

3.8 

.7 

17 

1.8  1.9 

•U 

ID 

2.6 

.5 

H* 

3.U 

.6 

18 

2.5  1.6 

.3 

U 

2.2 

.u 

1U 

2.8 

•5 

18 

An  interesting  uso  this  drag  data  can  be  £ut  to  is  the  rather  indirect 
neaairesent  of  tho  base  pressure.  The  dreg  coefficient  is  usually  broken 
up  into  the  head  drag,  skin  friction  drag,  and  base  drag  coefficients. 

(3)  *»-**♦  *Br  ♦ 


1.  The  Rlasiua  flat  plate  values  with  the  Van  Driest  correction  to: 
compressible  flow  were  used  for  the  lair.ihar  flow  computation  and 
tha  recendJy  derived' ■■formula-. obtained  by  Van  Drirst  yjJ  , DO 
wav  emplcwt*  foi  the  vtrb'Cent  boundary  lay*r.  The  Reynolds  numters 
baiau  on  diameter,  were* 

h Rd  x ID*5 

1.3  5*90 

1.8  8.17 

2.5  11.3U 

2.  The  percent  difference  vae  computed  with  respect  to  the  total  drag 

coefficient. 
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Slnns  %he  head  drag  arises  from  compressible  fluid  presauro  difference 
on  the  head,  it  can  be  computed  exactly  by  means  o t characteristics. 
This  vas  done  by  the  ZNIAfl  and  thi  results  are  given  in  Table  III. 

The  corresponding  values  obtained  by  ^he  linearised  the  cry  are  also 
listed  for  comparison  purposes.  In  figure  8 the  pressure  ratio,  as 
obtained  from  the  BfTAC  1«  plct-ted  against  distance  along  the  ncdsi. 


?A2Lg  HZ 

Head  Drag  feefficients 


M 

Acaet  Theory 

linearised  Theory 

1.3 

.071 91 

.0790 

1.8 

.0603 

.0660 

2.5 

.0>U9 

.£>  555 

An  examlra  tion  of  Table  III,  the  turbulent  values  of  Table  II  and 
figures  9,  ID,  -1  shews  that  contributes  about  202  of  the  drag 

while  and  7^  each  contribute  about  l£f.  The  ratio  of  bare  pressure 

to  free  e*uva*  pressure  can  now  be  computed  br  means  of  the  relation 

u)  £$  - 1 - jaLr  >wj 

y • 1 .U>5* 

Charters  and  Turatsly,  (j!3  J , were  able  to  make  an  independent  de- 
tendLnatlon  by  means  of  measuring  the  wake  angle  and  aaauatng  a 
Prandtl-Meyer  flow  around  the  base  to  thin  wake  angle.  A comparison 
of  these  methods  is  provided  in  Table  17.  T ho  values  for  Mach  number 
1.3  are  omitted  as  the  £NIAC  computations  of  ratio  of  pressure  at  end 
of  cylinder  to  free  strea  ware  made  only  as  low  as  Mich  number  1*5« 


1.  This  is  an  extrapolated  value  rs  was  (.ompuW  only  foi* 
If  9 1.?* 


U 
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TABLE  TV 


The  doscr  i-ar.;*  between  tha  tan  methods  la  caua»<Lby  the  neasureaent 
arror  1.;  the  wake  angle.  Tha  estimated  error  in  ¥B  tna  vaa*  angle  is 

*0 

.OU  whll*  the  vital  drag  method  has  as  estimated  arror  of  .01. 

An  a rami  nation  of  Tablo  IT  Indicates  that  does  not  vary  with 
the  length  of  tha  cylinder  and  hence  hsr  u the  only  varying  dreg  com- 
ponent. This  conclusion  is  further  substantiated  by  cowparing  the 
change  in  with  length  with  d hsr  in  Tablo  T.  The  change  in  Kp 


was  found  by  fitting  it  to  a linear  function  of  length  by  least  sou  ares 
and  d*DST  vai  found  by  differentiating  the  Tan  Driest,  formula. 


* ZVST  x 103  Slope  of  Kp  vs  L line  x ID3 

ft  ' L - 

.;  !*.o8  lull  + .li8 

.8  3.57  3-7U  ♦ .37 

.5  3.00  3.13*-.>> 

“no  ‘lr  at  digit  in  the  sc  dal  nurtV.  r indicates  its  length. 


TA&E  V 
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It  la  interesting  to  note  that  this  variation  decreases  with 
increa  lag  Mach  number. 

(b)  Spin  Deceleration  Mciient  Coefficient 

Tu<.  spin  decoloration  amaul  cue!  flu  lent  way  tn  defined  by 


(5)  K, 


p?tL^r 


* x 


where 

la  compare  Li  of  aerodynamic  acuct  alone  tho  eynmatry 

axis 

*Ld 

v • — — la  spin  in  radians  par  caliber 

<0^  is  ccmpamnt  of  angular  velocity 

Although  it  is  poaaiblo  to  determine  this  coef'lcient  from  the  jawing 
action,  after  part  of  the  program  had  bean  find  it  was  fan:!  th*v  this 
determination  van  not  very  accurate.  It  was,  therefore,  decided  to 
plaoa  Fins  in  tho  bases  of  the  remaining  missiles  and  measure  the  roll- 
ing motion  directly*-.  Figures  6 and  ? show  models  with  and  without  pins 
respectively.  The  individnal  round  values  of  are  tabulated  Ji 

Appendix  C and  plotted  in  Figure  12. 

Since  the  spin  deceleration  moment  is  a pure  viscous  effect,  it 
seemed  probable  that  it  could  be  related  to  Charters  and  Kent 

QlQ  have  shown  that  for  a cylinder  K&  * lA  i-ryy  *dis  appeared 

to  be  a good  approximation  for  our  configuration.  .(The  l/U  appears 
because  the  diaioter  and  not  the  radius  is  used  as  characteristic 
length).  Kpgy  was  obtained  from  Table  17  and  l/li  Kjgy  is  plotted  xn 

Figure  12.  The  agreement  seems  to  be  quite  good.  Since  the  absents 
ere  rneill  and  the  surface  conditions  from  round  to  round  are  clearly 
not  identical,  the  experimental  scatter  is  not  unexpected.  The  usual 
values  of  the  ladividuai  st*fj.rtxcal  standard  errors  as  given  in 
Appendix  C oeem  to  support  t hi  a consideration. 


1.  In  Appendix  A,  tho  exact  relations  wiployed  are  listed. 
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L.  Figure  13  is  clotted  .gainst,  length.  The  slopes  of  least 
squares  fitted  lines  are  .onpared  with  l/U  3KWF  in  Table  V. 


USZ  7 
Slope 


1.3 

.00106  * .00006 

.00302 

1.8 

.00099  ♦ .0000U 

.00089 

2.5 

.00079  *♦  .OOOOlt 

.00075 

TRANSTBtSE  F0BC2  AKD  “3KTVT  OOTmCWTS 

Tram  Appendix  A we  taro  the  statement  of  the  Xelley-HcShane 
linearity  aasueptinn. 

If  F,  and  7.  are  the  eoeponants  of  the  tnuuraerao  eero^jmaalo 
force,  then 

(A2)  F2  ♦ IF-,  - pd2**  £(-K„  + ivip  X ♦ (vljy  ♦ UCjJpTJ 

where  v la  the  spin  In  radlane  par  caliber 
V • Vj  ♦ IXj  la  the  coaplex  yaw 
<V  ieejd 

H - —lii..  a — la  the  com  lax  angular  Telocity 
*1 

and  the  I 'a  are  the  ballistic  coefficients.  They  are  identified  In 

both  AppeMix  A and  the  Table  of  Symbols  and  Coefficients.  Similarly 
for  the  eti.xments  of  the  tranerarse  moment  there  results 

(Alt)  Mj  ♦ i*}  “ (*-\2  f(-*Sr  “ ^ X * (**H  ♦ lw*n)*»3 

By  aseveing  a reaeooaole  siee  tar  the  Magnus  cross  epin  coefficients 
Kjy  and  K^,  ft  can  be  shown  that  they  hsT#  little  effect  on  the  aotiat 

of  m body  of  resolution  [23  j . In  none  a t'  the  firings  sale  on  the  BHL 
Cy^rk  Ranges  baa  this  aseurptior.  eeamed  unreraouable  aud  hence  oaV  the 
remaining  six  coefficients  will  be  considered. 
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In  Appendix  5 an  aerodynamic  nomenclature  is  defined  and  their 
relation  with  the  ballistic  coefficients  stated.  The.  conversion 
equations  are: 


hm  1 

*F  mTS 


* " i \ 

S “ ‘IS  Cn_ 


■ * 

! i 


«3'1S  <0»  *V 
q a 


•*  \ * V 


Not*  that  the  last  pair  of  relations  are  written  with  an  equality  sign 
and  not  In  the  Mannar  they  appear  in  (blD).  This  is  (tie  to  the  fact 
that  in  this  report  Kg  and  are  Measured  solely  by  Means  of  their 

effect  on  H in  tho  hawsgeaeous  part  of  equation  (AS)  and  >.a  is  stated 
in  Appendix  B the  equality  aipi  is  proper* 

Throughout  this  section  ns*  will  be  Bade  of  the  symmetry  argument 
that  all  ~ -efficients  are  even  functions  of  the  Magnitude  of  yaw. 

It  is  assumed  Jhat  che  yawing  notion  of  each  round  nay  be  characterised 
oy  its  mean  squared  yaw*  y , and  that  the  coefficient  obtained  fraa 
each  firing  la  constant  fur  ihe  yaw  encountered  in  each  flight  and  May 

V . ..ul  m 4h  T . fwnil  f Vi. . . (Iiinanrfiiinr. 


be  associated  wi  th  6 - It  wav  found  that  a s: 

“5 

on  0 was  sufficient  to  describe  the  data  for 


lep  1s  linear  dependence 
6Z<  30  square  decrees. 


(a)  Boreal  Tores  and  Overturning  Moment  Coefficients 

The  overturning  moment  coefficient  can  be  obtained  accurately 
Cron  the  turning  ratea  of  the  two  eras  of  the  characteristic  epic,  olio 
yawing  motion  of  a spinning  missile.  Actually  small  corrections  are 
necessary  which  involve  the  damping  exponents  and  the  Magnus  force 
coefficient  ty,  end  although  these  are  usually  lass  than  one  percent 

all  of  the  data  in  this  report  contain  them.  These  corrections  are 
explicitly  stated  In  Appendix  A.  A further  correction  is  necessitated 
by  the  variation  in  center  of  aass  froa  round  to  round*  Since  this 
variation  is  less  than  .02  calibers  for  models  of  the  same  type  all 
the  uudals  . re  corrected  to  a standard  center  of  uass  location  for  each 
type  by  weans  of  an  approximate.  in  equation  (A5*7). 

These  final  results  vers  exauined  for  yawing  notions  of  different 
amplitudes.  It  was  found  that  the  seven  caliber  forward  ckb.,  models 
together  with  the  nine  caliber  forward  and  middle  c.m.  models  had 
definite  dependence  on  the  magnitude  of  /aw.  Values  for  a fissd  Mach 


V-.--  •***«-• 


a linear  function  of  mean  squared  yaws 
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Tha  resulting  values  of  K*  _ are  listed  In  Table  VI  together  with 

V ' 

their  standard  errors t 


tABLK  VI 


V 

am 

Square  radians 

Mach  »«o. 

7 F 

9 F 

9 K 

1.3 

-36  ♦ 7 

-85  * 16 

-53  ♦ 26 

1.8 

-23  i 3 

-16  ♦ 3 

-59  ♦ 3 

2.5 

-53  ♦ 7 

-2D  ♦ 20* 

-56  ♦ 16 

»q.  (7) 

—1*1. 

-83. 

-39. 

“5 

* This  Value  ’ * poorly  determined  since  the  total  variation  of  ft 
Is  from  1 to  5.7  square  degrees. 


The  values  provided  by  Table  VI  should  be  eocstdered  only  good 
qualitative  results  since  they  were  obtained  by  a rather  crude  technique. 
Since  the  sables  are  mtall  the  standard  errors  have  only  qualitative 
significance.  With  this  la  mind  we  see  that  the  effect  of  yaw  on  tha 
overturning  moment  baa  a et-a  bill  sing  Influence  and  that  the 
of  these  Indications  of  non-linearity  lies  In  a non-linear  pressure 

distribution  over  the  rear  or  the  models1.  This  distribution  probably 
starts  about  six  calibers  frss  the  nose  and  increases  in  its  non- 
linear character  with  increasing  length.  From  this  we  would  expect  the 
longer  missiles  tn  exhibit  non-linearity  and  that  thia  effect  will  become 
more  pronounced  as  the  .enter  of  mass  is  moved  forward  and  thereby 
accentuates  the  effect  of  force  on  the  missile's  rear.  This  prediction 
is  roughly  verified  by  Table  VI. 

H.  a.  Sally  £y>”]  has  rscsntly  developed  a simple  relation  for 
Kjj  ^ end  IL  g by  consider  in*  the  viscous  cross  flow.  Uls  results  can 
6 6 

be  written  in  the  following  fomi 

(7)  2 - l./l*  C„  L2  p3r  - 21,] 

0 o 

(8)  K.,  „ - 3/1*  C-  L2 

"ft"  3c 


1.  Tills  ‘■▼plana tion  wan  suggested  by  J.  1>.  Uicoiaidej. 
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where  r Is  dlstanco  to  cantor  of  mass  from  the  nos*  and  Cp  la  ths  drag 

o 

coefficient  for  an  infinite  cylinder  In  a cruse  flow.  For  laalcar 
boundary  layers  Cp  is  1.2  and  for  turbulent  boundary  layers  iC  is  .25. 

Since  the  awran*  transition  point  is  about  1.5  cal.  nu-  of  Uia  ehuuldor 
far  tbs  seven  and  nine  caliber  rodels  a weighted  velhr’for  Ofj.  of  .78 

c 

was  need  for  the  seven  caliber  models  and  .66  for  the  nine  caliber  models. 
Kg  j nod  K^  2 have  been  computed  by  equations  (7)  - (6),  and  are  tabulated 
• y ■ g 

In  Tables  VI  and  Till  respectively.  Considering  tae  roughness  of  the 
experimentally  determined  ulues  and  the  use  of  a •weighted*  Cp  the 

*c 

agreement  is  good. 

-9 

Using  Table  VI  and  the  round  values  for  Kg  and  6 , individual 

values  of  L.  were  then  ooagmted  and  in  Figures  ?1*-16  are  plotted 
a 

against  Kacb  number.  The  marlimm  scatter  is  about  2%.  The  circled 
points  were  computed  by  fitting  £.  for  each  Mach  number  to  a line  as 

.o 

a function  of  eenter  of  aase  location:  (See  Sq.  A5-7).  The  slopes  of 

these  lines  are  the  normal  force  coefficients  and  are  tabulated  in 
Table  VII. 

TABLE  VII 

Normal  Three  Coefficient,  K^  ■ - -g  C^ 


M 

5 

7 

9 

1.3 

.98  ♦ 

.01 

1.02  ♦ ,02 

1.06  ♦ .oi 

1.8 

1.13  * 

.01 

i.13  ♦ .01 

1.16  ♦ .01 

2.5 

1.26  t 

.02 

1.21  ♦ .01 

1.28  ♦ .01 

Aw  is  dsacribed  in  Appendix  A it  la  also  possible  to  measure  Kg  by 

means  of  the  ;iwervinf  action,  of  were  obtained  by  this  means 

from  all  rounds  with  large  enough  swerving  motion2.  These  data  also 
e'  lbited  a dependence  on  yaw  for  the  seven  and  nine  caliber  rounds.  The 
values  of  K„  _ are  givv  n in  Table  VIII. 

V 

T.  More  precisely  /or  the  seven  caliber  models  Cp  • ^ ^.3/?a«. 

A sinilar  computation  was  perfomoc,  for  the  m&e  caliber  models. 

2.  Tha  criterion  was  that  tho  swerving  motion  be  six  times  the  linear 
measurenvirb  acruracy  of  .010-' . 
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sverr*  E^'a  vara  computed  in  a mannar  similar  to  thw  K^'a  and  ara 


plotted  ir.  Figure  17  • 


TABLE  mi 


V 

i 

Square  radians 

Mach  lo. 

7 

9 

1.3 

23  ♦ 7 

20  ♦ 3 
“ * 1 

1.0 

13  ♦ 3 

26  ♦ 3 

2.5 

l»6  ♦ 3 

33  ♦ 1? 

lq.  («) 

29. 

1*1. 

The  oirolau  value*  ara  those  obtained  froai  tha  Em' a at  different 

"o 

cantor  of  mass  locations  and  croaaaa  ara  baaad  on  toato  la  "tTL*  a 
Supersonic  Wind  Tunnel.  Tha  a greraant  la  fairly  good.  In  figure  18 
tha  diatanea  to  tha  orator  of  proMura  from  tha  noaa  in  calibara  C-  , 

rH 

la  plotted  a i obtained  frra  tha  individual  averring  ration,  tha  oantar 
of  aaaa  rathod,  and  tha  wind  Lunnol.  Sinca  this  acala  la  larger  than 
radal  a lac,  in  order  to  show  individual  points,  the  eurvea  are  r*- 
p lotted  againat  radal  length  at  the  actual  raoel  rlaa  of  20mn  par 
caliber  in  Figure  19.  Finally  Figure  30  ah  ova  X^'a  dapandenca  on 

length.  Hot*  that  both  cantor  of  proasuro  and  normal  force  are  ra- 
le tire  ly  insensitive  to  length.  (A  wild  exception  to  thia  la  the 
beharlar  of  the  cantor  of  preaeure  at  Kach  number  of  1.3)  In  thia 
characteristic  they  follow  tha  slander  body  prediction  that  cylindrical 
afterbodies  hare  no  affect  on  normal  fare#  or  cantor  of  pressure!. 

(b)  Magnus  Force  and  Moment  Coefficients 

In  thia  section  and  tha  next  oca  ve  will  discus*  thoaa  coefficients 
which  at  the  prasent  time  can  only  be  da  tormina  d by  the  precision  Range 
Technique.  Tha  determination  of  tha  quantities  la  difficult,  and, 
quit*  naturally,  la  not  as  accurate  as  Eg,  or  K,,.  The  difficulty 

lias  in  tha  fact  that  tha  dvnasdc  coefficients  effect  only  the  leaping 
of  the  epicyclic  yawing  action.  The  damping  is  difficult  to  determine 
a*  it  ia  toe  rate  of  change  of  a small  quantity.  Far  30mi  models  a 

1.  'The  data,  however,  is  not  very  close  to  tn*  slender  bo<fy  value*  of 
Ey  • .75  and  • i.iu.  This  discrepancy  has  been  observed  by 

many  investigator*  and  is  prooabiy  due  to  the  inappropriatenass  of 
the  theory. 
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reasonable  estimate  for  the  alia  of  arror  of  4 dirying  coefficient 

(L  la  2 x ID*5  1 /calibers.  Since  Xbr  the  fir*  call  bar  aodals 
* r 

25  x 10  , this  la  roughly  as  91  arror.  (Tha  figure  2 x w7 

la  sensitive  to  tha  freT-cr-7  of  tlx  t-^scds  % thalr  azplltud*^ 
and  tl'a  diatrtbution  rf  tha  eb.. 'rratian.) 

An  1 nap  action  -of  tha  Magrtua  boh  ant  ooaf  f*  slant  lata  for  tha 
TT,  ft,  and  9H  aodela  shows  that,  la  com on  with  tha  nofal  foroa 
and  overturning  eoeeni  eaafflolsnta,  thajr  in  inactions  or  yww.  H tla 
dapaodaiioo  la  oloodad  b y tha  aeattar  of  tha  data  and  uanea  K_  _ ai 

tabula  tad  in  Tabla  ZZ  la  poorly  dstandnad  for  soom  Mach  aaabara  and 
can  only  ba  astlaatad  for  othara. 

TABU  a 


V 

1 

Square  radians 

Maoh  Mg. 

TT 

ft 

m 

1.3 

-16  ♦ ID 

-1*3* 

- 16* 

1.8 

■43  ♦ 20 

-tt  ♦ 16 

- 1C" 

2.S 

-72  ♦ 2C 

•tt* 

-*• 

• Xrtimated  value 

*•  Obtained  trm  two  point* 

For  fixed  Mach  nuaber,  L la  k linear  function  of  0 an  tar  of  suae 
location  (saa  (A5.8)).  11  plotted  against  c.n*  In  Figures  21-23  and 

0 

tie  rather  largo  scatter  wariflea  tha  prsd'etlon  that  dynamic  data  la 
certainly  not  of  aa  good  quality  aa  that  for  Kg.  In  Table  1,  the 

a Lopes  of  tha  least  aquaraa  fitted  lines,  which  are  tha  Magnus  foroa 
joefflelenta,  are  tabulated.  Ibis  linear  relationship  la  used  to 
'btaln  'average*  nines  of  If  which  are  adhibited  aa  a function  of 

o 

Mach  rusher  in  Fi~ra  21*.  For  acne  rounds  vhsra  tha  yawing  action  la 
anM  that  '}••;  Mavnur  fores  has  a M.>aurabla  effect  on  tan  a. ere  lug 
notion  an  Independent  rceasarenant  of  Xj  la  possible.^  In  Figure  2$, 

Zy'  a aa  obtained  fras  tha  center  of  aais  eethod  tigetber  with  those 

obtained  free*  the  averring  notion  of  individuxl  rounds  ■ m plotted 
against  Haeh  numbers.  The  agreement  of  such  delicate  eessurenents  It 
quite  encouraging. 

1.  Appendix  A provides  the  -*xict  relations  eaplcyed  in  thia  datonainatlon 
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Hegnua  Force  Coefficient,  Kp,  From  C.H.  Method 


Mach  number 

1.3 

1.8 

2.5 

5 cal. 

.11,  ♦ .Cl 

.11.  * .0I4 

.05  ♦ .02 

7 cal. 

.20  ♦ .02 

.20  ♦ .02 

.US*  .02 

9 cal. 

.31  ♦ .03 

.33  ♦ .02 

. 26  ♦ .02 

«J.  C.  Martin  has  recently  suggested  a simple  nodal  by  which  Kp  and 

may  he  computed  theoretically  £2?^j  . Ha  shows  that  the  effect  of 

spin  on  a body  o i revolution  at  an  angle  of  attack  Is  to  rotate  the 
plana  of  symmetry  of  the  boundary  layer  configuration  slightly  out  of 
the  plane  of  yaw.  The  linearised  flow  oyer  the  resulting  shape  then 
provides  a force  which  is  perpendicular  to  the  plane  of  yaw  and  pro- 
portional to  the  magnitude  of  spin.  These  considerations  result  in 
the  following  formula  for  the  radio  of  slander  body  Magnus  force  coef- 
ficient Kp  to  slander  body  normal  force  coefficient  Kg  for  incompressible 

flow. 

<’>  £ . 6.91 't  «* 

where  L is  the  equivalent  cylinder  length  in  calibers1 

# 

and  A la  the  boundary  layer  displacement  thickness 
in  calibers. 


Since  our  main  interest  is  in  compressible  turbulent  boundary  layers, 
we  have  to  assume  that  Equation  (9)  will  apply  to  this  case.  From  [l3j 

and  Q3U3  » 11  oan  be  shown  that  a good  approximation  for  6*  is 

(l  - .U;  M)$  where  Rfl  is  th«  Reynold* ’ number 

baaed  on  the  diameter  and  c depends  on  the  velocity  profile.  Inserting 
w 

vhis  value  of  A together  with  the  oi. aider  body  normal  force  coefficient 


1.  Since  Martin's  calculation?  are  based  on  the  usual  flat  plate  boundary 
layer  assumptions,  his  vssults  are  for  a circular  cylinder  whose 
boundary  layer  build  up  ia  equivalent  to  the  body  of  revolution. 
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value  of  ^ in  (9)  'Jier.»  remuts: 


(ID) 


^ - c(l  - .11*  R.*1/5 


from  ( 10 ) it  la  nov  possible  to  compete  the  center  of  pressure  of 
Vie  Kagaua  Tore*  by  toe  formulas 

£ 


' ’TJ.  r 


„ *r 

. 9 £ 

Di 

where  £p.f  ia  the  equiva?.ent  cylinder  dj  stance  to  Ragnua  force  center 
of  pressure. 


We  now  assume  that  1),  the  equivalent  cylinder  length,  can  be  written 

L - L where  L le  the  modal  length  and  L is  a correct  wu  due  to  the  nose 
n n 

and  le  a function  of  Mach  number. 


(U*> 

(Ub) 


/• 

C*P.y  * C.P«F  ■ 1#^ 
d 

C.P._  • 9 I*  ♦ 5 !•„ 

F TZ  Tfi  a 

t - t(l  - .14  M)  (L  - L )9''  H.-l/5 


In  Figure  19  C.P..  is  plotted  against  length  and  it  can  be  seen 
r o 

that  the  predicted  elope  of  is  very  good.  For  fixed  Mach  number  wo 
fit  Equation  (11a)  and  obtain  the  following  values  for  hQt 


Mach  number  1.3  1.8  2.5 

ta  .52  ♦ .15  1.3*4  ♦ .08  2.55  * .12 

Using  equation  (lib)  as  an  empirical  relation  and  fitting  the  data 
ehown  in  lipuu  2r  there  results  c n .19  ♦ .03.  Those  requite  are  plotted 
in  Figure  17.  Thu  agreement  is  fair  for  ’such  a rough  theory,  although 
it  is  certainly  not  as  good  as  the  C.P.y  slope.  Finally  it  is  laportant 

note  tnit  the  slope  of  the  C.P.j  curve  versus  length  of  9/1 i*  is  in- 
dependent of  the  Mach  nunber  and  velocity  profile  and  depends  only  tn* 
the  assumption  that  for  t'irbulent  bcund.iry  laysrr  *7d  varies  inversely 

as  the  1/5  th  power  of  the  Reynolds  number  based  on  length. 


21 

RESTRICTCD— Se<y.-!ly  Information 


•wrvH" 


RESTRICTED—  Secwrdy 


Cc)  Damping  Force  and  Moment  Coefficients 


When  the  data  on  the  damptrg  »onant  coefficient  Kjj  ie  examined, 

once  egain  the  existence  of  non-linearity  in  the  7F,  9F,  and  SM  models 
can  be  seen.  In  Table  XI  i*  tabulated  the  Talus s of  which  were 

•uployed 

to  obtain  the  'a. 

TABLE  XI 

V 

Square  radians 

Mach  No. 

TF 

9F 

9H 

1.3 

590* 

650' 

650* 

1.8 

590  ♦ 361 

1*90  * 230 

720** 

2.5 

590  ♦ 230 

650* 

650* 

* estimated  **  twu  walnes  of 


Tnt”  ralies  are  Inserted  in  a sodlfied  Torn  of  (A5>9)  *nd  Kg,  the 
/«— fore*  coefficient  at  tb*  centroid,  is  obtained.1  The  da  aired 
fona  of  (X5>- 9)  ir 

[*h]  S < - <*M  " *H  * **3 

where  q is  measured  fro*  the  centroid  and  unatarred  quantities 
are  for  c.m.  at  centroid. 

In  Flgurea  26-28,  ^K^  ie  plotted  againet  c.*.  location  and  lines  art 
fitted.  The  ecatt»r  indicates  the  poor  quality  of  Kg.  In  order  to 
obtain  other  Taluea  of  Kg  for  c.».»s  which  are  not  at  tae  centroid, 
relation  (i5.lt)  sust  be  used.  Table  HI  present*  Kg  at  the  centroid 
and  it  is.  plotted  againat  Mach  number  in  figure  2?. 


1.  The  riddle  cents-  of  :'«se  rounds  hare  their  centers  of  ub  located  et 
their  geometric  centroid. 
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TABLE  XTI 
Kg  ct  Centroid 


Mach  Ho. 

5 

? 

9 

1.3 

-3.8  ♦ .1 

-5.3  * .6 

-8.f  * .8 

1.8 

-3.1  ♦ .5 

-5.1*  ♦ .1 

-9.2  *1,2 

2.5 

-1.1  ♦ .2 

-3.8  ♦ .1 

-6.6  *1.1* 

Iq.  (13) 

-U.6 

-7.0 

-9.3 

In  Figure  30  Kg  ax.  the  centroid  is  plotted  against  Mach  nuabar  while 
figures  31  and  32  give  the  centroid  values  of  Kg  and  as  functions 
of  length. 

In  order  to  get  soar  theoretical  basis  far  predicting  Kg  it  is 

necessary  to  use  certain  results  recently  obtained  by  Dor  ranee  fjSj  • 
According  to  Darrance's  "saro  order"  slender  body  theory  for  missiles 
without  boai.taii^.l  [The  relations  for  t4  and  C,,  were  first  obtained 
by  K.  Hnxk.] 

Sis  - (r  • '•* 


Ck  [tt  - r)!  • (r  - r.f] 


c„  - - 


■ 

U (L  - r) 


where  L is  length  in  calibers 

r is  distance  to  the  c.n.  from  the  nose  in  calibers 
v is  vo T nan  is  cal.'' 

r_  is  distance  to  the  centroid  from  the  ncse  in  calibers 
From  these  equations  there  results 


^ « (T.  - r)  - T 


T:  Care  has  to  bemused  in  order  to  transform  the  symbols  of  Q5J  to  t&ose 
of  this  report.  C and  C4  or  05^  actually  corrs spondee  C^  and 

nw  1 w & 

ai  d Cw  while  C and  C of  fjfl  correspond  to  C„  ♦ CM  and  CL.  * C,.  of  this 
m n«  m*  J »•.  w*  fi  ,i 

<r  a a 2j  q £ q a 
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Equation  (12)  predicts  tint  Kg  at  the  centroid  la  a pure  quadratic 


XW.cui.on  of  tiie  dioiai.ct.  ui  U.«a  centroid  from  the  base  which  has 


.79 


« It  a coefficient.  If  we  fit  Kj.  at  the  centroid  to  such  a function, 

there  results  the  values  of  1.19,  1.U3,  1.50  for  Kach  numbers  1.3,  1.8* 
and  2.5  respectively.  Since  the  fit  is  quite  good,  ve  will  consider 

Kg  • C (L  ro)2  as  a good  aaplrica?  formula  for  the  darp-ng  cnrer.t  at 

the  centroid.  Table  XII  compares  Kg  at  ce.  troid  with  Dorrance'a  Kg. 


The  agreement  is  not  too  satisfactory.  It  is  finally  of  intaxeat  to 
note  that  within  the  accuracy  of  this  elaeentary  theory  Kg  has  no 


contribution  froa 


idien  the  c.a.  is  at  the  centroid. 


DZMAKIC  STUBHITT 

I 

A study  of  the  dynamic  Stability  of  tha  rounds  fired  In  this 
program  gravidas  some  of  tha  most  interesting  results  of  this  report. 

As  in  [28]  a missile  is  defined  to  he  dynamically  stable  if  the  yawing 
notion  described  tqr  the  solution  of  the  homogeneous  saiation  of  the 
yawing  notion  does  not  in ere a ea.  It  is  proven  in  £?8j  that  a sufficient 
condition  for  $ namic  stability  of  a statically  unstable!  missile 
traveling  over  a flat  trajectory  is 

<1W  *-S-*D*‘&-k?,CA>0 

(u»b)  %7 - I 0 £ 1 K 2 

s (2  - s) 

A O ^ 

• . * (Gyroscopic  stability  factor) 

UB 

A " axial  moment  cf  inertia 
B " transverse  moment  of  inertia 
p ■ air  density 
d " disaster 

v • *“1^  ; spin  in  radians  per  calibe. 

“1 


1.  A missile  is  statically  unstable  if  0 . 
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(Dynamic  Stability  Factor) 

X 

A 

k^2  . (^  is  axial  radius  of  gyration  in  calibers) 

■ • mass 

2 

k£2  * (kj  is  transrsrse  radius  of  gyrction  in  calibers) 

jjf  j ^ q or  "i  y 2,  it  is  further  shown  that  a statically  unstable  ml  a si  la 

can  tot  be  stabilised  by  spin.  The  cures  • • 1 ■ is  plotted  in  Figure  33 

s(2  - s) 

and  the  stable  and  unstable  regions  are  identified. 


2(Kn  - Kd  - 

s — 1 — rj 

S - *d  ♦ k2\  - V 


This  reqairMant  ia  euch  sore  complete  than  the  classical  gyroscopic 
stability  requirement  that  |1L  Since  h ia  usually  posit  ire  for  missiles 

in  supersonic  flight1,  conditions  (lU)  will  reduce  to  the  classical  in- 


e quality  for  a ■ 1.  a is  tabulated  in  Appendix  C for  .11  roundo  and  it 
can  be  see  • the  e that  it  definitely  da  parte  tram,  this  optimum  raise  of 
unity.  For  soma  of  the  9?  modsla  a exceeds  two  and  hence  theaa  modala 
are  dynamically  unstable  aad  cannot  be  atabilisad  by  spinl  H baa  the 
further  property  that  the  slower  ana  has  the  stellar  daupinj^rate 
(c  ^ o^)  when  a 4 1 and  the  ravers#  la  true  (a 2*  o^,  whan  a>l. 


This  therefore.  Mans  that  the  ftster  arm  will  grow  for  the  9F  modala 
irrespective  of  apin.  (Unless  h ia  nagatira  the  Blower  arm,  however, 
will  always  shrink. ) 


It  la  possible  to  ****  an  important  jeensraliaation  of  the  dyii*»iu 
u -ability  in  tha  following  way,  (sea  Q4J  )• 


Theorem 


Tbs  damping  exponents  and  of  tha  apicyclio  yawing  artica  of  a 

statically  unstable  missile  are  greater  than  or  eqi'-al  to  an  assigned 
value,  « if  the  following  relations  are  satisfied: 

h - «>  0 

s * -±-= 

s(o)  [ Z - a(a)J 


1.  It  -as  positive  for  all  nissilas  llred  in  this  program. 
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C < ■ (oX  2 

whan  s(a)  - 2 (*N  ~ KP  m ’*1  *9  ~ J 

S - *P  ♦ - * 

(Mou  that  a(O)  ia  th<i  dynamic  Stufcili ty  factor.) 

With  tha  above  information  in  mind  wa  wall  new  move  on  to  a con- 
sideration of  tha  varloua  applications  o f our  experimental  information 
to  tha  stability  pro  bias.  First  tha  a f fact  of  cantar  of  mass  location 
la  of  interest* 

Looking  otar  tha  s' a for  tha  firs  caliber  models  in  Table  C-Zb, 
it  can  be  seen  that  *i  lias  between  .UD  and  1*30  and  hence  the  dynamic 
stability  is  of  little  intarsst  fbr  these  models.  The  dynamic  stability 
of  the  eaten  calibers,  sod  especially  the  nine  calibers,  is  more 
interesting.  In  order  to  obtain  a rough  picture  of  tha  effect  of  c.m. 
we  will  assume  that  c.m.  may  be  changed  while  the  naesse  and  radii  of 
gyration  remain  constant.  An  a rumination  shows  that  with  the  exception 
of  tha  bimetal  middle  c.m.'s  (9  M2,  9 Mi)  end  two  of  tea  rear  c.m.'s 
(9  HZ,  9 B3),  the?*  assumptions  are  roughly  true. 

With  these  assumptions  in  mind,  v2  required  for  stability  ia 
plotted  in  Tieure  3lt. against  c.m.  position  for  M ■ 1*3.  » was 
selected  as  one  *r  the  variables  in  tluj  plot  eo  that  tha  gyroscopic 
Stability  curve,  a • 1,  appears  as  a straight  line.  According  to 
this  plot  the  interval  of  c.m.  location  where  spin  stabilisation  ia 
possible,  identified  in  this  figure  fcy  • • Q.  ia  relatively  smell  end 
even  there  a rather  high  twist  is  required?  _The  rear  asymptote  c^-responda 
to  a - 0 while  tha  forward  one  is  caused  by  s - 2.  Mote  that_0.3  cal. 
rear  of  tha  centroid  is  located  at  tha  "optimum  point"  id»are  s • 1. 

Now  it  is  shown  in  £28]  that  tha  situation  is  improved  by  in- 
creasing fcj2.  Sines  k^2  ia  UPt  larger  for  the  bimetals  (9  M2  and  9 M3), 

tha  stability  curves  are  replotted  for  their  masses  and  radii  of  gyration 
in  Figure  35  (The  physical  constants  for  both  figures  ears  listed  in 
Table  mi). 


i 

■ ■ t 

1.  In  the  two  stability  plots  c.m.'s  forward  of  tha  centroid  are 
plotted  positively. 

i.  This  fact  was  fipat  observed  by  R.  Turetsky  in  fl8j  vhiph  was  an 
interin  report  on  the  program.  w 
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TABU  XIII 

m gins 

h2 

04 

•of* 

Standard 

137 

8.3 

.19 

51motal 

200 

7.6 

<26 

Tha  biaatala  as  can  be  easily  seen  art  much  more  stable  than  tha  solid 
models.  These  figures  become,  of  course,  more  inaccurate  aa  tha  c.a. 
is  mowed  fro*  tha  centroid.  If  tha  exact  physical  oharactaristlos  of 
tha  9 F's  are  usad,  a such  battar  determination  of  tha  c.a.  position 
for  which  s ■ 2 nay  ba  sad a*  Similarly  battar  Talus r fbr  tha  "optimum 
point*  and  tha  point  at  which  s ■ 0 can  ba  found  whan  tha  physical 
charactarisU.es  cf  tha  9 H'a  and  9 E*s  respectiwaly  are  used,  these 
calculations  have  bean  made  and  tha  results  appear  in  Table  HV.  Iota 
that  according  to  Table  XTV  it  ia  impossible  to  spin  stabilise  the 
9 F*s  and  fur  them  ora  that  at  M » 1.8  they  shoo  Id  ba  markedly  unstable. 
This  la  Tarifled  by  Table  C-Ub.  As  tha  yaw  increases,  howewer,  tha 
non-linearities  which  bare  bean  obserrad  throughout  the  prograa  han 
a destabilising  affect  and  increase  tha  else  of  a considerably* 

TABU  XT? 


Location  of  C.M.  Fro*  Centroid  For  Dynamically 
Stabla  Hina  Caliber  Models 


M 

s • 0 

■ ■ 1 (standard) 

a - 1 (bimetal) 

a • 2 

1.3 

-1.09 

-.32 

.05 

.97 

1.8 

-1.59 

-.70 

-.28 

.30 

2*5 

-2.03 

-.73 

-.28 

•9U 

As  a final  application  of  the  data  obtt  d by  this  prograa  we  will 
make  rr-ugh  estimates  of  the  stability  of  aodala  which  are  longer  thaw 
nina  calibers.  (Sines  these  estimate*  are  based  on  linear  theory,  they 
ara  at  ths  asrey  of  non-linearities  which  saw  to  increase  with  length.) 
For  simplicity  the  center  of  mass  will  be  located  at  ths  centroid.  Tha 
data  for  nine -caliber  length  module  definitely  shows  a bimetal  design  to 
be  superior  and  we  will  thus  consider  models  possessing  cylindrical  canter 
sections  o?  1»M*2th  2 J2„  and  of  density  different  from  tha  remainder  of 
the  modal.  Finally  the  following  formulas  for  centroid,  mass,  and 
moments  of  inertia  will  be  needed,  They  are  good  approximate  ens  for 
models  oxer  eight  calloer  lougi 
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c.m.  iron  bane  ■ l/'-'  (P  . ♦.?  ) 

, ■ length  ol  evlinder  in  calibers 
gi.  ■ if  (volume  of  nose  in  calibers)  • .8? 7 

■-  7*  (V3)  [Gh*  Gc4  &-1>21b] 

where  p1  is  density  of  nose  and  tail  material 
p2  is  density  of  center  section  material 

p-  $'+{*''**  2f»] 

where  • axLal  moment  of  inertial  of  ogival  nose 

■ .^565  (p^  d^) 

-fp  *d 

In  order  to  get  sonu  idea  of  the  stability  situation  for  these  long 
missiles,  we  will  select  fz  • 3 *nd  J?-  so  that  the  spin  required  'or 

Pi 

gyroscopic  stability  will  be  a mii.inur,  This  is  equivalent  to  requir.  lg 
that  be  a maximum.  In  addition  we  will  also  specify  p^  be  ecf\»l 
to  the  density  of  dural,  in  Table  XV  are  tabulated  the  resulting^,  A 
m,  and  s*s. 

TiBIA  XV 


T 

L (c'1>; 

.(.TS-c::!'") 

l.  (.r".-C?  L ) 

“ 

? 

2.U2* 

38 

3060 

31u 

1.0 

w 

3*l:5 

ns 

?080 

Jf-7 

1.2 

13 

3 . 

57 

36jO 

Uo5 

2.7 

16 

h,56 

72 

70o2 

55? 

3.2 
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According  t"  Table  IV,  models  with  c.a.  at  centroia  and  longer 
than  *' iven  calibers  *r?  lnpoaeible  to  atabilire  by  8pm.  It  Is 
possible  to  oaks  a rough  calculation  »a  to  tha  location  of  tha 
optioua  points  for  tba  models  In  Table  XV.  Thia  would  give  mm 
indication  of  possibla  improvement  of  dynamic  stability  by  varying 
c.a.  location  and  la  dona  in  ?abl«  XVI.  In  this  table  the  location 
of  tha  optima  point  ir  given  together  with  tha  gun  twist,  1_  , 

n 

required  for  tability. 

TABU  XVI 

I c.a.  rtiift  fraa  centroid  l/n 


16  -1.U6  1/9 

Tram  this  table  era  observed  the  very  important  facta  that  dynaalo 
stability  can  be  improved  ty  Moving  c.a.  rearward  and  that  'niita  long 
models  can  be  stabilised  with  reasonable  gun  twists,  me  rather  crude 
approximation  on  vhich  the  above  is  baaed  should  be  reempfcariaed  and  it 
should  be  remembered  that  the  above  is  done  only  ••  a rather  week  aid 
to  designers  of  longer  missiles. 


C.  ?/*• 

C.  R.  HORPHT  ' V 

£.  C.  J 


L.  X.  SCHMIDT 
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TA8L*  Of  SYMBOLS  COOTTCIBITS 


1 


*K 

b 

C 


(rXfr2,r3) 

o 

B 

* 

•i 

w 

*D8 

Ssr 

**V 

**T 


Adal  nouant  of  inartia 

Axial  kb*"1  flf  inartia  cf  tha  noaa 

frauaaaraa  aowat  of  l«ar-i* 

2 

a^lfloal  conataxit  daflnad  tflr  Bjj  6(X»  - Y0) 
Tho  dJmanaional  cylinder  drag  ooaffleiaat 

Magnus  foroa  cantar  of  proaiura 
Aorael  foroa  oanta.'  of  pws*»rn 

JA 

Aarodynaatc  foroa 

t'  * V * 1 * J T 

J*  “ *JD  ♦ k2t2’^ 

“®T"  11 

Spin  deceleration  aoaant  coaffiolant 

Drag  coaffiolant 

Axial  drag  ooaffioiaot 

Baaa  drag  coaffiolant 

Kaad  drag  eoeffioiant 

Skin  friction  drag  coaffiolant 

Mafsuj  f m*r«  coefficient 

Slr"d*»r  body  Magnus  foroa  coefficient 


1,  Hi  a symbols  which  appear  only  in  Ar-P*ndlx  3 ara 
table. 


omitted  from  this 
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VCMq  4 

N 

h 

VW 

VSle* 


Damping  moment  coefficient 


Lift  'arc#  coefficient 


Orcrt.urn.lng  (restoring)  moment  coefficient 


Noras!  fore*  coefficient 


MCN*  * 

MV 

pa 


Slander  body  normal  force  coefficient 
Dating  force  coefficient 


Magnua  moment  coefficient 


Magnus  cross  force  coefficient 
Magnus  cross  moment  coefficient 
Zero-yew  coefficient 
1 ^efficient  of  yew  scyiared  term 
Complex  constants  in  yew  equation 
Length  of  projectile 

Length  of  nose  which  is  unsffectire  in  Max-tin's  theory 

Equivalent  cylinder  length  in  Martin's  theory 
Mach  number 


Aerodynamic  moment 


Number  of  yaw  stations  and  timing  stations 


0 - 0* 


Shift  in  c.nu 
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fl-  V 


Const,*nt  term  in  tbs  mrra  equatl on 
Coefficient  of  linear  tarn  *a  the  averre  actuation 
Tho  "$*- function 

Tha  Reynolds  nuabar  based  on  disaster 


Tha  eoaffieienta  of  tha  exponential  terse  in  tha 
•xarra  aquation 

Constanta  dafinad  by  R^  ■ Ril  * *12 
Contribution  of  tha  lift  lores  to  tha  swart* 
Contribution  of  tha  Magma  fore*  to  tha  mm 

JL  “ ^2jT 

Tha  constant  tarn  in  tha  function 

Tha  coefficients  in  tha  polynomial  aanraasion  for  t 

Tha  coafflclant  of  - Jp  in  Rjj 

Tha  coefficient  of  tha  linear  tar*  in  tha  "Q=  function 
Tha  coefficients  in  tha  polynoadal  expression  for  0 


lha  coafflclant  of  J_  in  R. . 

3 9 _ 1 

Baplrical  constant  dafinad  by  L • c(l  - .U*)^  1 

* t * 

Tha  coafflclant  of  in  R^j 

Cantar  of  aasa  in  calibers  froa  tha  noaa 
Diameter 

Ths  coefficient  of  Jjj  in  R^ 

Acceleration  due  tu  gravity 
Vector  accalerutlon  ^ » to  gn-rLty 
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*? 

J?B 

u 

in 


h “ lD  ♦ ^2kH  ' h\ 

m ed2  ^ ig  *.hn  axial  radiua  of  gyration  in  calibers 

ed2  v 4,  1 1«  transverse  radius  of  gyration  in  cal:  oars 

" T ^ 

Ore-half  length  of  center  cylinder  in  bimetal  deu-igp 
Length  of  eyl'nder 

• ii  (velWof  nose)  Equivalent  length  of  nose 
a 


n 

n 


s 

s 

t 

u - (u1,u2>a3' 

v 

X 

y 

a 


Mass 

The  reciprocal  of  the  gun  twist 
Independent  variable 

Base  pressure 
Free  street*  pressure 

0-0*  (distance  flroe  centroid  in  definition  of  Sy* 

3i stance  in  calibers  fro*  the  nose  to  the  ce-oer 
of  mass 

Distance  in  calibers  fro®  the  noss  to  the  centroid 

Stability  factor 
Dynamic  stability  factor 
Tins 

Velocity  vector 

Vol'w.e 

Lateral  displacement 
Vertical  displacement 
Assigned  lower  limit  of  e^ 
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x-& 


♦ la.) 


(*2  ♦ i «j)d 


< 

(“1*  V ®3} 

( )• 

3!W 

~TT — 


r*w  draping  rates 

Ratio  of  specific  heat  of  air  rt  constant  prcoture 
to  that  at  constant  toIujm 

fa  - *i  v 4 - «i  sO  ^*1 

Magnitude  of  yaw  angle 
Mean  squared  yaw 

Boundary  layer  dieplacaient  thickness 

Perturbation  tens  in  solution  of  yaw  differential 
equation 

Standard  error  in  ballietic  coefficient 
Boll  angle 

Cowplex  yaw 

law  of  repose 

Complex  angular  velocity 
Spin  in  radians  per  caliber 

4 * 

Density  of  air 
Density  of  the  waterial 

Taw  turning  rates 

Angular  Telocity  of  the  missile 

Prints  denote  differentiation  vltn  respect  to  p 


Partial  deri’-ativ..  of  skin  friction  drag  coefficient 
with  respect  to  length 

Bars  denote  absolute  value 
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APMWDIX  At  SUMMARY  OF  THEORETICAL  RELATIONS  AH)  SFAT.K  RANQE 
TECHNIQUE 

In  the  bibliography  the  clasrical  publications  are  iioted  and  a 
fairly  complete  list  of  BR1>  publications  which  would  haws  a bearing  on 
this  ’■•port  is  provided*  Most  of  tha  theoretical  relations  used  hare 
•re  taken  fro*  [2k 7 or  _?J  * The  data  reduction  technique  la  described 

in  £lC]  or  in  the  Kopal  report  aentioned  i+\  the  introduction.  In 
this  appendix  ve  will  state  but  not  prove  the  relations  referred  to  in 

the  body  of  the  report.3  Although  these  relations  will  be  un  tenet  of 
tha  ballistic  K's  they  oan  be  easily  converted  to  the  aerodynamic  C'e 
where  necessary  by  uae  of  Appendix  B. 

We  first  da  fin a a right  handed  orthogonal  coordinate  37a  tea  with 
axes  numbered  1,  2,  3 ecring  with  the  eiaalla  and  so  orientated  that  tha 
1 axis  always  points  along  tha  missile's  axis  and  tha  2 axis  lias  in 
tha  b orison tal  plana  and  points  to  tha  right,  lbs  linear  velocity  of 
the  center  of  use  and  angular  velocity  of  the  missile  are  sspnesaed 
in  this  coordinate  system  aa  tha  vectors  (u^,  Ug,  u^)  and  (m^,  »2,  a^) 

respectively.  Iwploylng  the  convenient  represents ticn  of  cceqtlax 
variables  the  basic  XellayHteBhane  linear  force  system  far  a missile 
possessing  an  ar  jle  of  rotational  symmetry  lass  than  120°  and  a plans 
of  mirror  eyraetry  is  defined  by  the  following  equations: 

(Al)  rx  - ^ { ft  - - pd2  **  Xp 

(A2)  T2  * ir3  - pd*u*  £(-*„♦  >■  ♦ Mjy  ♦ l*s)J 

(A3)  Kj^  - -pd3  u 12  vEa 

(Ah)  Kj  ♦ iMj  • pd3u,2  H-vEy  - U^)  1 ♦ (-5g  ♦ ivEjj)!*] 


(F, , Fj,  T^i  la  the  aerodynamic  force  vector 

F is  n-wpanent  of  aerodynanic  force  vector  directed  along  the 
trajectory 

(li^,  My,  Kj)  is  the  serr dynamic  moment  vector 
p is  the  density  of  air 
d Is  the  diameter  of  the  node! 

1.  Although  all  of  the  relations  stated  are  trua  for  all  missiles 
possessing  the  proper  symmetry  required  for  (Al)  * (AU),  soma  of 
the  remarks  ou  ranpe  technique  refsr  only  to  spinning  bodies  of 
revolution.  30 
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is  the  oucnpiex  yaw 


to,d 

v • ~-t—  is  nondimenoional  spin 
(m,  ♦ iax)d 

H ^ — i«  complex  angular  teloeity 

Kjj  ftrag  Fares  Coefficient1 

’Su  Axial  Drag  Fores  Coefficient 

Spin  Decslsration  Moment  Coefficient 
Kjj  Normal  Force  Coefficient 

£y  Magnus  Force  Coefficient 

Kjp  cross  Spin  Mag.au  Force  Coefficient 
*3  Damping  Foret  Coefficient 

Kj,  Magma  Moment  Coefficient 

Kjj  Overturning  Moment  Coefficient 

Iq  Damping  Moment  Coefficient 

Kgp  Cross  Spin  Magnus  Moment  Coefficient 

Since  the  cements  are  defined  with  respect  to  center  of  mass  and 
X is  defined  from  the  motion  of  the  center  of  mass,  then#  ballistic  K's 
are  functions  cf  center  of  mass  location.  Since  vs  require  symmetric 


1.  For  xiro  cross  spih  it.  car  be  ea«ily  shown  that  Kg  ’•*  cee  g ♦ 

i KDA  + ^ “ Kp4  wb«,r*  6 " * X | 

* 


K^O  sin  6 
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mass  distribution,  *11  c.a. 1 s are  on  the  axis.  T»s  explicit  form  of 
these  .unctions  1st 


(A?.l)  - kd 

UfJ.2)  K„#  - K,, 

U5.3)  1/  * *F 

(AS.U)  Eg*  - Eg  ♦ 

(AS.*)  Ijf  - Kj,  ♦ 


(45-6)  Ka*  - Ka 

(AS.7)  V * *M  '*3 
(4$. 6}  tj*  - tj  * O, 

(45.9)  «a*  - «h  ♦ q(lg  ♦ V4  *\ 
(45.1&)K1I*  - Ij,  ♦ q(Kjy*  Ej)*  q^ 


♦^•rs  ths  starred  quantities  correspond  to  tbs  center  mss  located  at  point 
O*  and  tbe  una tarred  quantities  to  esntar  of  nasa  at  0.  q la  the  axial 
distance  frost  0 to  0*  in  calibers^  and  ia  considered  positive  when  Measured 
toward  the  beet  fro*  the  nose. 

Placing  these  definitions  in  the  equations  of  notion  we  obtain  the 
following  equations  for  ths  axial  and  yawing  notion  for  a flat  trajectory? 

4 


(46) 

(47) 


"Jr 


Dv 


(48)  X»  ♦ (H  - iT)  X*  ♦ (-*  - iwT)  X - 0 

where  primes  denote  differentiation  with  respect  to  non  dimensional 

y* 


axial  arc  length  p 


■r 


and  t is  tine. 


JD-t 


m.  m j* 

kf*  ■ —j*  (k^  is  axial  radius  of  ^ration  in  calibei  s) 

4 id  axial  meaent  of  inertia 
n is  the  mass 


1.  4 caliber  is  a unit  of  length  equal  to  the  nissils' a dianwter. 

2.  Equation  (Afl)  is  baled  on  the  sire  assumption  or  convention 
that  « tsims  may  be  omittsd  in  compariaor  vlth  J teims. 

I 

uO 
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-2  »d2 

■ -g-  (kj  la  tranavsrao  radius  of  gyration  in  calibers) 
B is  Jim  transverse  moment  of  inertia 


" • Vs  (J»  ♦ »2  S*  V * ‘2'!  * 

* - •’»  - - V*  'r 

0 • T'  - O'.*  k?"2  V * SJ  T 

C*2  * «i*p  4 l^«3  * «lVJ  d 

f 5 

“i 

(gj,  g2,  gy  la  t>M  vector  accalaratloa  da*  to  gravity 
*2  ♦ Oj  - X 

Tba  solution  to  the  aquation  of  yawing  notion  can  bs  wrlttsn  in  tbs 
form 

(»)  1 . ^ .‘--*1  * *1*  . I,  t *4*  . V, 

stare 

K'lsM 

Xj  art  caaplex  constants  da  ponding  on  initial  conditions 

V,  - - (yaw  of  repose) 

11  u^cM 

g la  acceleration  doa  *.  ■»  gravity 

are  constants  and  {^'  are  linsar  functions  of  p 
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Tha  ojcponente  of  (49)  nay  ba  ralatad  to  tha  coefficient*  of  equation 
(48)  u/i 

(410)  v - ♦ tf2’ 

(411)  U - ♦ «2>  - *t\ 

» H'  •*>*  -v»*  £0 

T--V, 

D - 


U12) 

(413) 

(Alii) 


V ♦ V 


»l'  * *2 


whsre  a very  good  appro  trijaatlaa  to  tha  parturbatlon  tara(f^  la 


mi 


- A alanlla  It  aald  to  ba  atatioalljr  atabla  if  0.  A statically 

yft 

unstable  adsalle  la  aald  to  ba  eyroacoplcally  atabla  If  a - f l.  A 

adaaila  la  dynamically  atabla  If  tha  ya*  daacrlbad  by  tha  hoaofaaaoo a 
aolutlon  to  (48)  doaa  not  lnaraaaa*  Tor  a statically  tmatabla  aiaalla 
tha  aaponanta  a^  and  a^  ara  graatar  than  or  equal  to  an  aaaignad  tralna 

a If 


(41$)  H ♦ D - a > 0; 


(416)  efi 


and 


■(*  - *) 

(417)  0 < a < 2 

whara  a (a)  - -jj? c la  tha  ganarallaad  dynamic  stability  factor. 

If  o • 0,  a become  tha  dynamic  stability  factor  a (0)  and  (41$)  - (41^) 
bscoan  conditions  for  dynaale  stability.  If  a (0)  doaa  not  aatlafy 
(417),  a statically  unstable  aiaalla  can  never  ba  dynamic ally 

atabiliaed  by  spin1. 


1.  For  a statically  stable  mioeila  (41$)  and  (416)  apply  vhan  a la 
outslda  tha  interval  TO,  H abila  only  (41$)  is  needed  when 
a is  inside. 
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1C  the  position  of  tha  missila  is  calculated  from  tha  aquations  of 
motion,  and  tha  p axis  is  taken  to  ba  down  range  in  tbs  horirontal 
plana,  tha  y axis  pointing  up,  and  tha  x axis  detcr»io»d  by  tha  right 
hand  rule,  we  have  tha  relation 

(A10)  x*ly»?  + Q?*  s(hi1  4 1(V)p  ♦ RjKg  s(“*2  4 

*ffh^rJT>  * u dp  dp 


where 


P and  Q are  ccaplex  constants  da  tern  load  by  Initial  conditionfl 
x and  y are  in  ealibera 

Ri  “ *11  * **12 

Ril"  ‘il  * bU  JT  * °L1  JS  * dil  JXT 


«V2  - «i  ) 


2v  ^ 04* 


' (J?  r?>*  (vv  * w77?>’  * 


J.  ♦ 


V*,' 


" 5 ■■■  "p  gr» 

* «i  > 

Ri2  * *12* JN  ’ JJp  * bi2  4 ci2  J5  * *12  JXT 


#l'  ®i  s *(^'2  - •t2)  «> 


♦ & 
(iv7  ♦ «!*) 


JS  ‘ 


Jxy 


Tha  first  integral  is  the  displacement  due  to  yaw  of  repose  and 
can  usually  ba  estimated  to  a sufficient  accuracy  for  range  work,  lha 
real  part  of  this  integral  ia  callsd  the  •drift'1.  Tha  second  integral 
is,  of  coui'se,  the  gravity  drop.  The  expression  can  na  replac'd 

by  J^,  the  lift  force  coefficient. 


i 
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Or  the  ERL  Spark  Photo graohy  Rang*  tha  drag  coefficient  la  faund 
by  star,  of  up  to  twaiva  tir-e-diatance  aaamiraeents.  Tha  dlataoca 
arror  la  leas  than  >001  feet  and  tha  laaat  count  for  tha  tine  la  5/8 
of  a nicrosecond.  Tha  data  ara  usually  flttad  to  a cubic  in  diatanca. 


(A19)  t • aQ  ♦ a,p  ♦ a^p2  ♦ a^p^ 


where  t ie  tine 


<A*>)  i 


Tha  velocity  u at  point  * la  than  glean  by 
1 


♦ 2 ajp  ♦ 3 a^p 


and  JQ  at  point  p can  than  bu  cospuUd  fron  (l6). 
2al  + 604) 

(A21)  Jp  - r 


♦ 2^  ♦ J*jP 


« and  Jp  ara  usually  evaluated  at  tha  cantar  of  tha  data.  Tha  toraparatura 

and  prassura  are  naasurad  bafora  aaeh  firing,  tha r* by  ; -oviding  tha 
Telocity  of  sound  *nd  dansity  of  air.  Ihon  this  tha  Mach  minbar  and  tha 

ocP  nay  ba  ccnputad.  (a.  A,  B,  and  ran  tar  of  aaaa  location 


dansity  factor 


ara  precisely  Masurod  for  aaafa  nodal  bafora  firing). 


0 can  ba  diraetly  datanainad  by  naaauring  tha  spatial  location  at 
aach  atation  of  too  pint  placed  in  tha  baaa  of  tha  nodal.  This  t.'n.i 
datarainea  tha  roll  angle  9 aa  a function  of  position,  p,  on  tha  rang*. 
Ttaaaa  data  ara  thau  flttad  to  a cubic  polysonlal 


(A22)  9 ■ b0  ♦ bj^  p ♦ bjp2  ♦ b^p^ 


From  (A7)  we  haver 

2bj  ♦ 6bjP 


(A23)  D - 


bj.  ♦ 2b  jP  ♦ 3*^, 


Proa  (A23)  and  (A 21)  can  than  obtain  J^. 


Tbn  two  conponanta  of  tha  yaw  of  the  nijsile  ara  usually  meaaured 
to  an  accuracy  of  .001  radians.  They  are  tlien  fitted  by  a conbinatlaa 
graphical  and  analytical  technique  to  equation  (A?)*  Froa  tha  eoef* 
ficiarta  af  thia  fit  by  naana  of  (All)  - (All*)  wa  can  ootain  H,  T, 

anJ  a rvlatiuely  poor  aecond  determination  of  D.  Tha  spatial  position 
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of  tha  jenter  of  mass  ic  now  measured  to  an  accuracy  of  .001  fact  and 
fitted  to  aquation  (118).  Of  tha  fern  umrally,  however,  only 

ia  Mall  determined.  Fbrtunatoly  - Jp  is  the  principal  constituent 

for  moat  firing a and  c-n  than  be  determined.  Thia  plua  tha  valued  of 
og.  H,  a.JI  T tl.ri.  proy.ua  ua  with  Jg  and  In  certain  cares  Eg?  n*n 

be  data  rained,  and  from  thia  follow  valuta  of  Jy.  Gy  firinr  different 

canter  of  mass  positions,  equations  (AS)  than  provide  ns  with  Js»  V 

and  a ascend  determination  of  Jg. 

In  au—ry  we  are  that  the  firing  of  a single  nodal  with  satis- 
factory initial  yawing  notion  (large  enough  to  measure  and  null  enough 
td  be  linear)  and  satisfactory  averring  motion*  will  provide  values 
Ijj,  Xg,  Xjj,  Kj,  Xg,  and  possibly  Xy  at  a given  Mach  number,  firings 

of  identical  models  with  different  center  of  mass  positions  at  tha  law 
Mach  cumber  then  yield  additional  values  of  Kg,  Kj,  end  Xg. 


1. 


It  can  easily  be  shown  that  the  lift  force  contribution  to  tha  slow 
Nerve  arm  is  (Rj)l  " JXj  U *21  * *22  JL  and  the  careen. on^pf  con- 
tribution tram  the  Magnus  force  is  (l^y  ■ | Kg  * bjj,  y. 

for  most  rounds  J3  and  Jjry  may  be  osd.ttadja|.><,  Sj^nd  bj-^  *y/J-  a^^  J. 


Hence  (Rg)^  '—JXpJ  j RjjJ-  Swerving  motion  is  satisfactory  for 
Kg  when  (Rg)L  is  more  than  aim  tines  the  experimental  accuracy . Since 
cur  Magnus  force  measures  were  limited  in  number,  Kj,  raluos  are  con- 
sidered when  (Rj)p  is  twice  the  experimental  accuracy. 
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APPhTTOIX  Bi  CONVERSION  OF  THE  BALLISTIC  COEFFICIENTS  TO  AERODYNAMIC 

coefficients 

The  work  of  Able  report  has  been  dor\e  in  toms  of  tho  ballistic 
K'«  which  oro  little  known  outside  tho  field  of  bolllotics  and  ujr  bo 
quits  confusing  to  ah  aerodynwiciot  who  does  his  dynamic  stability 
analy*ea  in  tens  of  tho  aerotynAmlc  C'a.  Zt  is  thorofors  worthwhile 
V express  tho  rj«ult*  of  this  report  in  toms  of  those  symbols.  This 
effort  is  handicapped  however,  by  tho  three  faetsi 

1.  The  missiles  usually  treated  in  ballistics  have  a rotation 
symmetry  which  results  in  pairs  of  aerodynamic  coefflolente 
being  equal  and  hence  corresponding  to  only  one  ballistic 
coefficient. 

2.  In  ballistics  tho  missiles  usually  have  a high  rate  of  spin 
and  Magrrj.  affects  have  to  be  considered  to  which  there  are 
no  corresponding  aerodynamic  coefficients. 

3.  Terns  involving  tbs  rate  of  change  of  angle  of  attack  appear 
in  most  aerodynamic  stability  analysis  while  no  such  terms 
appear  lu  the  usual  ballistic  force  system. 


The  swlal  components  of  the  aerodynamic  force  and 
usually  defined  In  aerodynamic  nomenclature  aai 


int  are 


X • - 1/2  P V2  S C„ 

L • 1/2  P T2  5 b($5>  Cj 


where  p la  air  dere^ty 

V is  axial  velocity 
S le  a reference  area 
b is  the  wing  span 


Fran  this  we  see  that 


- 1/2  S/d  C, 


- - 1A  (SbVd1*)  C jf 

v - ^ (pb/2V) 
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xt  the  traiureraa  cop.por.sntt  of  the  aerodynamic  force  end  Mcaient 
art  assumed  to  be  linear  functions  of  yiw,  change  in  yaw,  and  angular 
Telocity,  and  Magnus  coupling  it  introduced,  wa  have  tha  following 
definitional 


(B3) 


(1/2 

u V2  S)  j 

[%  ’ ' 

■ % <N>  ♦ 

% $D 

* J*r 

, * 4 % 

J $>] 

(V* 
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/s 

H <U>  • ' 
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:2,  <#] 

*rv 

1 ♦ C- 

V 
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<$>]  <$A 

(1/2 

p 8 c) 

I 
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* C«,  (5T*] 
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*[V 
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‘♦0 

<CT>  * cv 
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> <=.,  $} 

*[v 

» 
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••s. 
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Where  c it  tha  wing  chord  and  the  angles  a,  P and  angular 
valaoitias  A,  f,  q and  w art  those  defined  in  tha  standard  aero- 
dynamic ncsaanolatura. 


If  tha  Missile  is  aeauMad  to  oossass 
symaatry,  it  follows  froa  Qlt  J 

trigonal  cr  greater  rotational 

L2!/ 

% ‘ % * \ 

-C-  • c,  ■ c„ 

P*  pp  "pa 

c*  - c,  J 

”pq  pr  * pq 

% « “ % ’ \ 

Cr  “ Cz  .(|)S  Cj. 
'p*  " "pA 

U7 
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The  third  Mt  of  ayubola  ia  introduced  in  order  to  emphasise  the 
existence  of  aymmatry  end  will  be  employed  throughout  the  remainder 
of  this  appendix.  If  we  insert  these  symbols  into  (S3),  multiply  the 
second  and  fourth  equations  by  i and  add  to  the  first  and  third  respec- 
tively these  re  suit at 

I ♦ il  • (1/2  p V2  S)  (o.  ♦ 1 <R>  °v](p  * 1 *> 

' ‘E*  V 1 \] 

(B5)  H ♦ IS  • (1/2  p f « S)  f - ‘ °»J  <»  ♦» -> 

•[81  V - 

If  equation  (£6)  is  compared  with  aquations  (A2)  and  (Ah),  the  Magma 
end  tion -Magnus  static  coefficients  are  easily  related. 

- - 1/2  S/d2  Cj, 

ft 

- 1/2  Sc/d3  O, 

"a 

ty  - 1/U  Sb/d3  t 

" P° 

» 

Kj  » - 1/U  Scb/d4 
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lie  re iationWiips  between  the  remaining  tynaaic  coefficients  is 
somewhat  more  uur.plicatad.  Fortunately  it  can  easily  be  shown  that  the 
remaining  Ha  gnu  a coefficients  are  lost  in  the  differential  equations  of 
jrawln g Motion  due  to  tlie  <T  contention.  It  therefor#,  remains  only  to 
connect  two  ballistic  coefficient.*,  Sg  and  F^,  with  four  aerodynamic 


cosffj  clouts,  Cjj  , CN  , and  . 

q A q A 


>.d*: 


Tn  order  to  do  this  we  need  only  to  consider  the  purpose  of  this 
work,  namely  to  state  the  —suits  of  this  report  in  aerodynaaic  nomen- 
clature. Sines  this  report  Is  concerned  with  stability,  the  only  con- 
tribution of  the  aerodynamic  coefficients  is  how  they  appear  in  the 
basic  differential  aquations.  This  Mans  that  in  order  to  obtain  ths 
partner  of  Ig  we  see  what  coefficient  appears  in  ths  corresponding 

point  of  the  different!#!  equation  similar  to  (l)  which  Is  based  on 
the  aerodynamic  force  system  (see  f 271  for  example).  Gy  thia  tactic 
we  hares  U J 


(37). 


c2S 


\ * V 


Since  the  -raj  or  function  of  Kq  Is  its  contribution  to 
center  of  mass  is  altered  we  base: 


when  tbs 


(c, 


♦ cM  ) 
q *4 


Notei  The  method  of  obtaining  (G7)  and  (B6)  is  not  toe  de  air  able. 
It  would,  of  course,  be  mors  satisfying  to  enlarge  the  ballistic  fores 
system  so  that  there  would  exist  a one-to-one  correspondence.  It  also 
should  be  noted  that  (87)  follows  from  a comparison  of  the  homogeneous 
equations. g In  the  yaw  of  -epose,  equation  (A9),  should  be  replaced 

toy  - lA  -r~  Cj,  . 
d4  "q 


By  use' of  (32),  (B6),  (B7)  end  (Bfl)  it  is  now  pose*  ole  to  conrert. 
our  symbols.  We  will  merely  tabulate  the  result*.  (t^  will  be  replaced 

bry  Kjj  - Kp  for  thia  purpose.) 


'I  i 


- [c»a  * 2 % * v*  v!  <!>*  * <*•>] 


*-3<R>4 


M 


H - k 


-2  ocS 


In  order  tc  avoid  contusion  c«  ini  (a*)  will  be  replaced  by  y** 

k9 
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;-w,‘  v^vv-^v7 

for  bodies  of  resolution  if  S is  the  iaxiao*  cross-sectional  ere* 

S * l*  *2  “ 5 * * * 1 

.\S“I  S>  % 

S"i  S#  J *nmi  \ 

(BID)  Ky  - 16  CNpa  5 S - ‘ 16 

S-*1?  (sa  4 v » s-*  -£  \ V ■ * 

.ayvK’sj** 

• ■ zffirzr 1 1 ’ v vK'vv- K" v T* 

Z~*Z  i-^,s7^  | * L wh-rs  L is  to.  aodal  length  in  calibers. 

This  je lection  seems  to  complicate  our  equations  uinscessr-ily. 
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APPENDIX  Cs  TABLES  OF  DATA 

In  Table  0-1,  the  physical  characteristics  of  the  twenty-seven 
no  del  types  are  tabulated.  The  types  are  Identified  by  two  numerals 
separated  by  a letter.  The  first  number  gives  the  model  length,  the 
letter  «pecifi*s  whether  the  center  of  mass  is  forward  (T),  Middle 
(M),  or  rear  (P.«,  ami  t v second  number  identifies  different  ly<«  ■? 
of  the  sane  length  and  canter  cf  mesa  location.  The  compost lien  of 
each  Model  is  given  by  three  letters  which  soecify  the  metals  used 
in  the  nose,  center,  and  base  sections  respectively. 

In  Tables  C-2,  C-3,  and  C-U,  thp  aerodynamic  date  for  each  round 
ere  given.  The  mean  squared  yaw,  6*,  is  in  square  degrees  and  is 
effectively  aero  for  those  rounds  where  it  is  omitted.  The  dreg 
coefficient,  X~,  ia  tabulated  for  all  rounds  possessing  over  five 
timing  stations,  and  the  spin  decelerating  coefficient,  K^,  for  only 

rounds  with  pins.  Values  of  Kg  are  given  for  those  rounds  for  which 

both  arms  of  their  epicyulic  yawing  motion  emceed  .00$  radians. 

Values  of  Kg  and  K»  are  listed  when  both  aims  exceed  .007  rauians. 

In  addition  to  the  arm  else  requirement  there  must  be  fifteen  observa- 
tions and  a fivc^sble  distribution  of  the  observations  on  the  epicycle. 
For  ease  reminds  it  was  possible  to  calculate  Kg  from  the  spin  and  the 

turning  vats  of  <--ae  arm  when  only  ona  arm  exceeded  .00$  radians  in 
aims  end  the  uodel  possessed  pins.  Kg  was  calculated  when  the  swerve 

associated  with  it,  was  greater  then  .06  inches.  For  those 

rounds  idiich  did  not  have  KD  or  Kg  values.  Kg  and  ij  were  computed 

using  values  corresponding  to  the  asms  type  at  the  some  Mach  nuab er. 

The  column  marked  N - Rj  gives  the  total  number  of  observations 
and  the  number  of  time  measurements,  and  are  In  l/caliber*  and 

v is  in  radiame/caliber.  v uay  be  converted  to  gun  twist  It  n by  the 
relation  n • yS  . 

In  Table  C-$,  the  aerodynamic  data  for  models  possessing  mean 
squared  yaws  of  over  thirty  square  degrees  can  be  found.  Table  C-6 
gives  values  of  Magnus  force  coefficients  Measured  from  the  swerving 
motion  of  those  models  who.«e  Magnus  sworving  notion,  (Rj)y,  ia  greater 

thwt  ,02  inches,  d.tno  tie  Jter.isticai  error  of  the  -arioue  ballistic 
^•v.fficients  was  fairly  uniform  for  nods!  types  possessing  the  same  spin, 
only  representative  values  are  given  in  Table  C-7*  Finally  the  tumiiu? 
rates  of  the  two  epicycle  an\.a  are  provided  by  Table  C-6. 
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In  '•.one lusion  the  numbering  syebom  for  the  ®^6J* 
described.  This  can  be  cone  oy  the  following  table  for  the  five 

caliber  length  model*. 


5*01  to  $-29 
5-31  to  5-59 
5-61  to  5-69 
5-91  to  5-99 


Forward  c.m. 
Midale  C.m. 
gear  c.». 
Large  yaw 


Tb.  numbers  for  the  7 and  9 caliber  length  are  dlrided  t»  the 
way. 
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